King Sejong’s National Referendum on Tax Reform (ggplot)

<<<<<<< HEAD ======= >>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

Data

원자료는 세종실록, 요약표는 오기수 교수의 논문에서.

지역	계급	관료여부	찬반	집계
서울	대신 등	관료	찬성	21
서울	대신 등	관료	반대	194
서울	3품이하현직	관료	찬성	259
서울	3품이하현직	관료	반대	393
서울	3품이하전직	관료	찬성	443
서울	3품이하전직	관료	반대	117
유후사	품관촌민	품관촌민	찬성	1123
유후사	품관촌민	품관촌민	반대	71
경기	수령	관료	찬성	29
경기	수령	관료	반대	5
경기	품관촌민	품관촌민	찬성	17076
경기	품관촌민	품관촌민	반대	236
평안	대신 등	관료	반대	1
평안	수령	관료	찬성	6
평안	수령	관료	반대	35
평안	품관촌민	품관촌민	찬성	1326
평안	품관촌민	품관촌민	반대	28474
황해	수령	관료	찬성	17
황해	수령	관료	반대	17
황해	품관촌민	품관촌민	찬성	4454
황해	품관촌민	품관촌민	반대	15601
충청	대신 등	관료	반대	2
충청	수령	관료	찬성	35
충청	수령	관료	반대	26
충청	품관촌민	품관촌민	찬성	6982
충청	품관촌민	품관촌민	반대	14013
강원	수령	관료	찬성	5
강원	수령	관료	반대	10
강원	품관촌민	품관촌민	찬성	939
강원	품관촌민	품관촌민	반대	6888
함길	대신 등	관료	반대	1
함길	수령	관료	찬성	3
함길	수령	관료	반대	14
함길	품관촌민	품관촌민	찬성	75
함길	품관촌민	품관촌민	반대	7387
경상	수령	관료	찬성	55
경상	수령	관료	반대	16
경상	품관촌민	품관촌민	찬성	36262
경상	품관촌민	품관촌민	반대	377
전라	대신 등	관료	반대	2
전라	수령	관료	찬성	42
전라	수령	관료	반대	12
전라	품관촌민	품관촌민	찬성	29505
전라	품관촌민	품관촌민	반대	257

전체 찬반

# pander(format(tbl.vote, big.mark = ","))
kable(t(as.matrix(format(tbl.vote, 
                         big.mark = ","))), 
      caption = "전체 찬반", align = 'c')

전체 찬반

찬성	반대
98,657	74,149

kable(t(as.table(format(tbl.vote, 
                        big.mark = ","))), 
      caption = "전체 찬반", 
      align = 'c')

전체 찬반

찬성	반대
98,657	74,149

kable(t(as.table(format(prop.table(tbl.vote)*100, 
                        digits = 3, 
                        nsmall = 1))), 
      caption = "전체 찬반(%)", 
      align = 'c')

전체 찬반(%)

찬성	반대
57.1	42.9

원형그래프

<<<<<<< HEAD

source("./theme_kr.R")
=======
source("./theme_Malgun_kr.R")
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
# str(theme.kr)
str(tbl.vote)
##  'xtabs' int [1:2(1d)] 98657 74149
##  - attr(*, "dimnames")=List of 1
##   ..$ vote: chr [1:2] "찬성" "반대"
##  - attr(*, "call")= language xtabs(formula = counts ~ vote, data = sejong.poll.kr)
kable(t(as.matrix(tbl.vote)))



찬성
반대




98657
74149



vote.df <- as.data.frame(tbl.vote)
n.vote <- length(levels(vote.df[, 1]))
y.coord <- vote.df$Freq/2 + c(0, cumsum(head(vote.df$Freq, -1)))
pie.label <- paste(levels(vote.df$vote), 
                   format(vote.df$Freq, 
                          big.mark = ","), 
                   sep = "\n") 
pie.label
## [1] "찬성\n98,657" "반대\n74,149"
kable(vote.df, 
      align = rep("c", 2))



vote
Freq




찬성
98657


반대
74149



str(vote.df)
## 'data.frame':    2 obs. of  2 variables:
##  $ vote: Factor w/ 2 levels "찬성","반대": 1 2
##  $ Freq: int  98657 74149

p1 <- ggplot(data = vote.df, 
             mapping = aes(x = "", 
                           y = Freq, 
                           fill = vote)) 
p2 <- p1 + 
  geom_bar(width = 1, 
           stat = "identity", 
           position = position_stack(reverse = TRUE))
p2
<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
p2.2 <- p1 + 
  geom_bar(width = 1, 
           stat = "identity", 
           position = "dodge")
p2.2
<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

p3.2 <- p2.2 + 
  theme_bw() + theme.kr
p3.2
<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
p4.2 <- p3.2 + 
  scale_x_discrete(name = "찬반") +
  scale_y_continuous(name = "집계", 
                     breaks = vote.df$Freq, 
                     labels = format(vote.df$Freq, big.mark = ","))
p4.2
<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

p5.2 <- p4.2 +
  scale_fill_manual(name = "찬반", 
                    values = rainbow(n.vote)[n.vote:1], 
                    labels = vote.df$vote)
p5.2
<<<<<<< HEAD

p6.2 <- p5.2 +
  geom_text(aes(y = vote.df$Freq/2), 
            label = format(rev(vote.df$Freq), big.mark = ","), 
#            family = "HCR Dotum LVT", 
            position = position_dodge(width = 1))
p6.2

=======

p6.2 <- p5.2 +
  geom_text(aes(y = vote.df$Freq/2), 
            label = format(rev(vote.df$Freq), big.mark = ","), 
#            family = "Malgun Gothic", 
            position = position_dodge(width = 1))
p6.2

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

p7.2 <- p6.2 +
  theme(legend.position = c(0.25, 0.75))
p7.2
<<<<<<< HEAD

p8.2 <- p5.2 +
  guides(fill = "none") +
  geom_text(aes(y = vote.df$Freq/2),
            family = "HCR Dotum LVT", 
            label = rev(pie.label),
            position = position_dodge(width = 1))
p8.2

=======

p8.2 <- p5.2 +
  guides(fill = "none") +
  geom_text(aes(y = vote.df$Freq/2),
            family = "Malgun Gothic", 
            label = rev(pie.label),
            position = position_dodge(width = 1))
p8.2

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
# ggsave("../pics/sejong_geom_bar_total_ggplot.png", p7.2, dpi = 72)

p3 <- p2 + 
  theme_bw() + theme.kr
p3 
<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
p4 <- p3 + 
  scale_x_discrete(name = "찬반") +
  scale_y_continuous(name = "집계", 
                     breaks = cumsum(vote.df$Freq), 
                     labels = format(cumsum(vote.df$Freq), big.mark = ","))
p4
<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

p5 <- p4 +
  scale_fill_manual(name = "찬반", values = rainbow(2)[2:1], labels = vote.df$vote)
p5
<<<<<<< HEAD

p6 <- p5 + 
  geom_text(aes(y = y.coord), 
            label = format(vote.df$Freq, big.mark = ","), 
#            family = "HCR Dotum LVT", 
            position = position_identity())  
p6

=======

p6 <- p5 + 
  geom_text(aes(y = y.coord), 
            label = format(vote.df$Freq, big.mark = ","), 
#            family = "Malgun Gothic", 
            position = position_identity())  
p6

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
# ggsave("../pics/sejong_geom_bar_total_ggplot_stack.png", p6, dpi = 72)

p7 <- p6 +
  theme(legend.position = c(0.2, 0.75))
p7
<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
p8 <- p5 + 
  guides(fill = "none") +
  geom_text(aes(y = y.coord),
            label = pie.label,
<<<<<<< HEAD
            family = "HCR Dotum LVT", 
            position = position_identity())
p8

=======
            family = "Malgun Gothic", 
            position = position_identity())
p8

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

pie.total.1 <- p2 + 
  coord_polar(theta = "y", 
              start = 3 * pi /2 , 
              direction = -1)
pie.total.1

<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

pie.total.2 <- pie.total.1 + 
  scale_y_continuous(name = "", 
                     breaks = NULL) +
  scale_x_discrete(name = "") 
pie.total.2

<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

pie.total.3 <- pie.total.2 +
  scale_fill_manual(name = "", 
                    values = rainbow(n.vote)[n.vote:1])
pie.total.3

<<<<<<< HEAD

pie.total.4 <- pie.total.3 +
  theme_void(base_family = "HCR Dotum LVT")
pie.total.4

=======

pie.total.4 <- pie.total.3 +
  theme_void(base_family = "Malgun Gothic")
pie.total.4

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

pie.total.5 <- pie.total.4 +
  guides(fill = "none")
pie.total.5

<<<<<<< HEAD

pie.total.6 <- pie.total.5 +
  geom_text(aes(y = y.coord), 
            label = pie.label, 
            family = "HCR Dotum LVT", 
            position = position_identity())
pie.total.6

pie.total.7 <- pie.total.6 +
  ggtitle("전체 찬반") +
  theme(plot.margin = unit(c(1, 1, 1.5, 1), "lines"), plot.title = element_text(hjust = 0.5, size = 20))
pie.total.7

ggsave("../pics/sejong_total_pie_ggplot.png", pie.total.7, dpi = 72)

=======

pie.total.6 <- pie.total.5 +
  geom_text(aes(y = y.coord), 
            label = pie.label, 
            family = "Malgun Gothic", 
            position = position_identity())
pie.total.6

pie.total.7 <- pie.total.6 +
  ggtitle("전체 찬반") +
  theme(plot.margin = unit(c(1, 1, 1.5, 1), "lines"), plot.title = element_text(hjust = 0.5))
pie.total.7

# ggsave("../pics/sejong_total_pie_ggplot.png", pie.total.7, dpi = 72)

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

pie.gg()

pie.gg <- function(df, ggtitle = "", font.family = ""){
  n <- length(names(df))
  y.coord <- df$Freq / 2 + c(0, cumsum(head(df$Freq, -1)))
  pie.label <- paste(levels(df$vote), format(df$Freq, big.mark = ","), 
                     sep = "\n") 
  p1 <- ggplot(df, aes(x = "", 
                       y = Freq, 
                       fill = vote)) 
  p2 <- p1 + 
    geom_bar(width = 1, 
             stat = "identity", 
             position = position_stack(reverse = TRUE))
  pie.1 <- p2 + 
    coord_polar(theta = "y", 
                start = 3 * pi / 2, 
                direction = -1)
  pie.2 <- pie.1 + 
    scale_y_continuous(name = "", 
                       breaks = NULL) +
    scale_x_discrete(name = "") 
  pie.3 <- pie.2 +
    scale_fill_manual(name = "", 
                      values = rainbow(n)[n:1])
  pie.4 <- pie.3 +
    theme_void(base_family = font.family)
  pie.5 <- pie.4 +
    guides(fill = "none")
  pie.6 <- pie.5 +
    geom_text(aes(y = y.coord), 
              label = pie.label, 
              family = font.family, 
              position = position_identity())
  pie.7 <- pie.6 +
    ggtitle(ggtitle) + 
    theme(plot.margin = unit(c(1, 1, 1.5, 1), "lines"), 
          plot.title = element_text(hjust = 0.5))
  return(pie.7)
}
dump("pie.gg", 
     file = "./pie.gg.R")

pie.gg(vote.df, 
       ggtitle = "전체 찬반", 
<<<<<<< HEAD
       font.family = "HCR Dotum LVT")

======= font.family = "Malgun Gothic")

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

계급 및 지역별 찬반

계급별 찬반

kable(format(vote.class, 
             big.mark = ","), 
      align = "r", 
      caption = "계급별 찬반")

계급별 찬반

	대신 등	3품이하현직	3품이하전직	수령	품관촌민
찬성	21	259	443	192	97,742
반대	200	393	117	135	73,304

kable(format(prop.table(vote.class, 
                        margin = 2)*100, 
             digits = 3, 
             nsmall = 1), 
      align = "r", 
      caption = "계급별 찬반(%)")

계급별 찬반(%)

	대신 등	3품이하현직	3품이하전직	수령	품관촌민
찬성	9.5	39.7	79.1	58.7	57.1
반대	90.5	60.3	20.9	41.3	42.9

품관촌민 별도 분석

품관촌민의 수효가 상대적으로 많아서

kable(format(vote.class.2, 
             big.mark = ","), 
      align = rep("r", 2), 
      caption = "품관촌민")

품관촌민

	관료	품관촌민
찬성	915	97,742
반대	845	73,304

소계를 교차표의 주변에 계산

kable(format(vote.class.2.am, 
             big.mark = ","), 
      caption = "관료와 품관촌민(소계)", 
      align = rep("r", 3))

관료와 품관촌민(소계)

	관료	품관촌민	계
찬성	915	97,742	98,657
반대	845	73,304	74,149
계	1,760	171,046	172,806

백분율을 계산하여 주변에 집계.

kable(format(prop.table(vote.class.2, 
                        margin = 2)*100, 
             digits = 3, nsmall = 1), 
      caption = "관료와 품관촌민(%)", 
      align = rep("r", 3))

관료와 품관촌민(%)

	관료	품관촌민
찬성	52.0	57.1
반대	48.0	42.9

data frame 변환

vote.class.2.df <- as.data.frame(vote.class.2)
v.names.class.kr <- c("찬반", "계급", "집계")
kable(format(vote.class.2.df, 
             big.mark = ","), 
      align = c('c', 'c', 'r'), 
      col.names = v.names.class.kr, 
      caption = "관료와 품관촌민")

관료와 품관촌민

찬반	계급	집계
찬성	관료	915
반대	관료	845
찬성	품관촌민	97,742
반대	품관촌민	73,304

vote.bureaus.df <- subset(vote.class.2.df, 
                          vote.class.2.df$class.2 == "관료", 
                          select = c("vote", "Freq"))
kable(format(vote.bureaus.df, 
             big.mark = ","), 
      align = 'r', 
      col.names = v.names.class.kr[-2], 
      caption = "관료의 찬반")

관료의 찬반

찬반	집계
찬성	915
반대	845

vote.commons.df <- subset(vote.class.2.df, 
                          vote.class.2.df$class.2 == "품관촌민", 
                          select = c("vote", "Freq"))
kable(format(vote.commons.df, 
             big.mark = ","), 
      align = 'r', 
      row.names = FALSE, 
      col.names = v.names.class.kr[-2], 
      caption = "품관촌민의 찬반")

품관촌민의 찬반

찬반	집계
찬성	97,742
반대	73,304

원형그래프

pie.bureaus <- pie.gg(vote.bureaus.df, 
                      ggtitle = "관료의 찬반", 
<<<<<<< HEAD
                      font.family = "HCR Dotum LVT")
pie.commons <- pie.gg(vote.commons.df, 
                      ggtitle = "품관촌민의 찬반", 
                      font.family = "HCR Dotum LVT")
=======
                      font.family = "Malgun Gothic")
pie.commons <- pie.gg(vote.commons.df, 
                      ggtitle = "품관촌민의 찬반", 
                      font.family = "Malgun Gothic")
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
pies.grid <- grid.arrange(pie.bureaus, 
                          pie.commons, 
                          ncol = 2, 
                          top = "")

<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

pies.grid

## TableGrob (2 x 2) "arrange": 3 grobs
##   z     cells    name                 grob
## 1 1 (2-2,1-1) arrange       gtable[layout]
## 2 2 (2-2,2-2) arrange       gtable[layout]
<<<<<<< HEAD
## 3 3 (1-1,1-2) arrange text[GRID.text.1265]

ggsave("../pics/sejong_bureaus_commons_pie_ggplot.png", pies.grid, width = 8, height = 4, dpi = 72)

======= ## 3 3 (1-1,1-2) arrange text[GRID.text.1227]

# ggsave("../pics/sejong_bureaus_commons_pie_ggplot.png", pies.grid, width = 8, height = 4, dpi = 72)

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

지역별 찬반

관료와 품관촌민

kable(vote.region.bureaus, 
      caption = "관료들의 지역별 찬반")

관료들의 지역별 찬반

	서울	경기	평안	황해	충청	강원	함길	경상	전라
찬성	723	29	6	17	35	5	3	55	42
반대	704	5	36	17	28	10	15	16	14

kable(format(prop.table(vote.region.bureaus, 
                        margin = 2) * 100, 
             digits = 3, 
             nsmall = 1), 
      align = rep("r", 9), 
      caption = "관료들의 지역별 찬반(%)")

관료들의 지역별 찬반(%)

	서울	경기	평안	황해	충청	강원	함길	경상	전라
찬성	50.7	85.3	14.3	50.0	55.6	33.3	16.7	77.5	75.0
반대	49.3	14.7	85.7	50.0	44.4	66.7	83.3	22.5	25.0

kable(format(vote.region.commons, 
             big.mark = ","), 
      align = rep("r", 9), 
      caption = "품관촌민들의 지역별 찬반")

품관촌민들의 지역별 찬반

	유후사	경기	평안	황해	충청	강원	함길	경상	전라
찬성	1,123	17,076	1,326	4,454	6,982	939	75	36,262	29,505
반대	71	236	28,474	15,601	14,013	6,888	7,387	377	257

kable(format(prop.table(vote.region.commons, 
                        margin = 2) * 100, 
             digits = 1, 
             nsmall = 1), 
      align = rep("r", 9), 
      caption = "품관촌민들의 지역별 찬반(%)")

품관촌민들의 지역별 찬반(%)

	유후사	경기	평안	황해	충청	강원	함길	경상	전라
찬성	94.1	98.6	4.4	22.2	33.3	12.0	1.0	99.0	99.1
반대	5.9	1.4	95.6	77.8	66.7	88.0	99.0	1.0	0.9

서울의 찬반

kable(vote.seoul.class, 
      caption = "서울의 찬반")

서울의 찬반

	대신 등	3품이하현직	3품이하전직
찬성	21	259	443
반대	194	393	117

kable(format(prop.table(vote.seoul.class, 
                        margin = 2) * 100, 
             digits = 1, 
             nsmall = 1), 
      align = rep("r", 3), 
      caption = "서울의 찬반(%)")

서울의 찬반(%)

	대신 등	3품이하현직	3품이하전직
찬성	9.8	39.7	79.1
반대	90.2	60.3	20.9

막대그래프

Stack

vote.seoul.df <- as.data.frame(vote.seoul.class)
n.fill <- length(levels(vote.seoul.df[, 1]))
x.stack <- vote.seoul.df[, 2]
y.stack <- unlist(tapply(vote.seoul.df$Freq, 
                         vote.seoul.df$class, 
                         function(x){x / 2 + 
                             c(0, cumsum(head(x, -1)))}))
y.breaks <- unlist(tapply(vote.seoul.df$Freq, 
                         vote.seoul.df$class, 
                         cumsum))
y.label <- format(y.breaks, big.mark = ",")
b1.seoul <- ggplot(vote.seoul.df, 
                   aes(x = x.stack, 
                       y = Freq, 
                       fill = vote)) +
  geom_bar(stat = "identity", 
           position = position_stack(reverse = TRUE))
b1.seoul

<<<<<<< HEAD

b2.seoul <- b1.seoul +
  theme_bw(base_family = "HCR Dotum LVT") +
=======

b2.seoul <- b1.seoul +
  theme_bw(base_family = "Malgun Gothic") +
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
  theme.kr +
  scale_x_discrete(name = "계급") +
  scale_y_continuous(name = "집계", 
                     breaks = y.breaks, 
                     labels = y.label) +
  scale_fill_manual(name = "찬반", 
                    values = rainbow(n.fill)[n.fill:1], 
                    labels = vote.seoul.df$vote)
b2.seoul  
<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

b3.seoul <- b2.seoul +
# geom_text(aes(y = y.stack/2), label = vote.seoul.df$Freq)
  geom_text(aes(y = y.stack), 
            label = vote.seoul.df$Freq, 
            position = "identity")
b3.seoul
<<<<<<< HEAD

ggsave("../pics/sejong_seoul_barplot_stack_ggplot.png", b3.seoul, width = 6, height = 4.5, dpi = 72)
=======

# ggsave("../pics/sejong_seoul_barplot_stack_ggplot.png", b3.seoul, width = 6, height = 4.5, dpi = 72)
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

Dodge

n.fill <- length(levels(vote.seoul.df[, 1]))
b1.seoul.dodge <- ggplot(vote.seoul.df, 
                         aes(x = x.stack, 
                             y = Freq, 
                             fill = vote)) +
  geom_bar(stat = "identity", position = "dodge")
b1.seoul.dodge

<<<<<<< HEAD

b2.seoul.dodge <- b1.seoul.dodge +
  theme_bw(base_family = "HCR Dotum LVT") +
=======


b2.seoul.dodge <- b1.seoul.dodge +
  theme_bw(base_family = "Malgun Gothic") +
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
  theme.kr +
  scale_x_discrete(name = "계급") +
  scale_y_continuous(name = "집계", 
                     breaks = vote.seoul.df$Freq, 
                     labels = vote.seoul.df$Freq) +
  scale_fill_manual(name = "찬반", 
                    values = rainbow(n.fill)[n.fill:1], 
                    labels = vote.seoul.df$vote)
b2.seoul.dodge  
<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

#  N <- nrow(vote.seoul.df)
#  index <- as.vector(matrix(1:N, nrow = 2)[2:1, ])
y.label <- unlist(tapply(vote.seoul.df$Freq, 
                                 vote.seoul.df$class,
                                 rev))
b3.seoul.dodge <- b2.seoul.dodge +
  geom_text(aes(y = vote.seoul.df$Freq/2), 
#            label = vote.seoul.df[index, "Freq"], 
            label = y.label,
            position = position_dodge(width = 0.9)) +
  ggtitle("서울의 찬반")
b3.seoul.dodge
<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
# ggsave("../pics/sejong_seoul_barplot_dodge_ggplot.png", b3.seoul.dodge, width = 6, height = 4.5, dpi = 72)

Fill

n.fill <- length(levels(vote.seoul.df[, 1]))
y.fill <- unlist(tapply(vote.seoul.df$Freq, 
                        x.stack, 
                        function(x){cumsum(x)/sum(x)}))
p.fill <- unlist(tapply(vote.seoul.df$Freq, 
                        x.stack, 
                        function(x){(x / 2 + c(0, cumsum(head(x, -1))))/sum(x)}))
b1.seoul.fill <- ggplot(vote.seoul.df, 
                        aes(x = x.stack, 
                            y = Freq, 
                            fill = vote)) +
  geom_bar(stat = "identity", 
           position = position_fill(reverse = TRUE))
b1.seoul.fill

<<<<<<< HEAD

b2.seoul.fill <- b1.seoul.fill +
  theme_bw(base_family = "HCR Dotum LVT") +
=======


b2.seoul.fill <- b1.seoul.fill +
  theme_bw(base_family = "Malgun Gothic") +
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
  theme.kr +
  scale_x_discrete(name = "계급") +
  scale_y_continuous(name = "집계(%)", 
                     breaks = y.fill, 
                     labels = format(y.fill * 100, 
                                     digits = 2, 
                                     nsmall = 1)) +
  scale_fill_manual(name = "찬반", 
                    values = rainbow(n.fill)[n.fill:1], 
                    labels = vote.seoul.df$vote)
b2.seoul.fill  
<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

b3.seoul.fill <- b2.seoul.fill +
  geom_text(aes(y = p.fill), 
            label = vote.seoul.df$Freq, 
            position = "identity") +
  ggtitle("서울의 찬반")
b3.seoul.fill
<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
# ggsave("../pics/sejong_seoul_barplot_fill_ggplot.png", b3.seoul.fill, width = 6, height = 4.5, dpi = 72)

barplot.gg

Stack

barplot.gg.stack <- function(df, 
                             base_family = "", 
                             ggtitle = "", 
                             xlab = "", 
                             ylab = "", 
                             fill.name = ""){
n.fill <- length(levels(df[, 1]))
x <- df[, 2]
y <- unlist(tapply(df$Freq, 
                   x, 
                   function(x){x / 2 + c(0, cumsum(head(x, -1)))}))
y.breaks <- unlist(tapply(df$Freq, 
                         x, 
                         cumsum))
y.label <- format(y.breaks, big.mark = ",")
<<<<<<< HEAD
delta <- (max(y.breaks) - min(y.breaks)) / 20
y.breaks.sort <- sort(y.breaks)
diff(y.breaks.sort) < delta 
index <- which(diff(y.breaks.sort)  > delta)
y.breaks <- c(0, y.breaks.sort[c(index, length(y.breaks.sort))])
=======
# delta <- (max(y.breaks) - min(y.breaks)) / 20
# y.breaks.sort <- sort(y.breaks)
# diff(y.breaks.sort) < delta 
# index <- which(diff(y.breaks.sort)  > delta)
# y.breaks <- c(0, y.breaks.sort[c(index, length(y.breaks.sort))])
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
y.label <- format(y.breaks, big.mark = ",")
b1 <- ggplot(df, 
             aes(x = x, 
                 y = Freq, 
                 fill = vote)) +
  geom_bar(stat = "identity", 
           position = position_stack(reverse = TRUE))
b2 <- b1 +
  theme_bw(base_family = base_family) +
#  theme.kr +
  scale_x_discrete(name = xlab) +
  scale_y_continuous(name = ylab, 
                     breaks = y.breaks,
                     labels = y.label) +
  scale_fill_manual(name = fill.name, 
                    values = rainbow(n.fill)[n.fill:1], 
                    labels = df$vote, 
                    guide = guide_legend())
b3 <- b2 +
  geom_text(aes(y = y), 
            label = format(df$Freq, 
                           big.mark = ","), 
            position = "identity") +
  ggtitle(ggtitle)
return(b3)
}

Dodge

barplot.gg.dodge <- function(df, 
                             base_family = "", 
                             ggtitle = "", 
                             xlab = "", 
                             ylab = "", 
                             fill.name = ""){
n.fill <- length(levels(df[, 1]))
x <- df[, 2]
y.dodge <- df$Freq
b1 <- ggplot(df, 
             aes(x = x, 
                 y = Freq, 
                 fill = vote)) +
  geom_bar(stat = "identity", 
           position = "dodge")
b2 <- b1 +
  theme_bw(base_family = base_family) +
#  theme.kr +
  scale_x_discrete(name = xlab) +
  scale_y_continuous(name = ylab, 
                     breaks = y.dodge,
                     labels = format(y.dodge, big.mark = ",")) +
  scale_fill_manual(name = fill.name, 
                    values = rainbow(n.fill)[n.fill:1], 
                    labels = df$vote)
#  N <- nrow(df)
#  index <- as.vector(matrix(1:N, nrow = 2)[2:1, ])
y.label <- unlist(tapply(df$Freq, 
                         df[, 2],
                         rev))
b3 <- b2 +
  geom_text(aes(y = y.dodge/2), 
#            label = format(df[index, "Freq"], big.mark = ","), 
            label = y.label,
            position = position_dodge(width = 0.9)) +
  ggtitle(ggtitle)
return(b3)
}

Fill

barplot.gg.fill <- function(df, 
                            base_family = "", 
                            ggtitle = "", 
                            xlab = "", 
                            ylab = "", 
                            fill.name = ""){
n.fill <- length(levels(df[, 1]))
x <- df[, 2]
y.fill <- unlist(tapply(df$Freq, 
                   x, 
                   function(x){cumsum(x)/sum(x)}))
p.fill <- unlist(tapply(df$Freq, 
                        x, 
                        function(x){(x / 2 + c(0, cumsum(head(x, -1))))/sum(x)}))
b1 <- ggplot(df, 
             aes(x = x, 
                 y = Freq, 
                 fill = vote)) +
  geom_bar(stat = "identity", 
           position = position_fill(reverse = TRUE))
b2 <- b1 +
  theme_bw(base_family = base_family) +
#  theme.kr +
  scale_x_discrete(name = xlab) +
  scale_y_continuous(name = ylab, 
                     breaks = y.fill,
                     labels = format(y.fill * 100,
                                     digits = 2,
                                     nsmall = 1)) +
  scale_fill_manual(name = fill.name, 
                    values = rainbow(n.fill)[n.fill:1], 
                    labels = df$vote, 
                    guide = guide_legend())
b3 <- b2 +
  geom_text(aes(y = p.fill), 
            label = format(df$Freq, big.mark = ","), 
            position = "identity") +
  ggtitle(ggtitle)
return(b3)
}

barplot.gg

barplot.gg <- function(x, 
                       position, 
                       base_family = "", 
                       ggtitle = "", 
                       xlab = "", 
                       ylab = "", 
                       fill.name = ""){
  switch(position,
         stack = barplot.gg.stack(x, 
                                  base_family = base_family, 
                                  ggtitle = ggtitle, 
                                  xlab = xlab, 
                                  ylab = ylab, 
                                  fill.name = fill.name),
         dodge = barplot.gg.dodge(x, 
                                  base_family = base_family, 
                                  ggtitle = ggtitle, 
                                  xlab = xlab, 
                                  ylab = ylab, 
                                  fill.name = fill.name),
         fill = barplot.gg.fill(x, 
                                base_family = base_family, 
                                ggtitle = ggtitle, 
                                xlab = xlab, 
                                ylab = ylab, 
                                fill.name = fill.name))
}

(bar_seoul_stack <- barplot.gg(vote.seoul.df, 
                               position = "stack", 
<<<<<<< HEAD
                               base_family = "HCR Dotum LVT", 
=======
                               base_family = "Malgun Gothic", 
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
                               ggtitle = "서울의 찬반", 
                               xlab = "계급", 
                               ylab = "집계", 
                               fill.name = "찬반"))

<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

# ggsave("../pics/sejong_seoul_barplot_stack_ggplotv2.png", bar_seoul_stack, width = 6, height = 4.5, dpi = 72)

(bar_seoul_dodge <- barplot.gg(vote.seoul.df, 
                               position = "dodge", 
<<<<<<< HEAD
                               base_family = "HCR Dotum LVT", 
=======
                               base_family = "Malgun Gothic", 
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
                               ggtitle = "서울의 찬반", 
                               xlab = "계급", 
                               ylab = "집계", 
                               fill.name = "찬반"))

<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

# ggsave("../pics/sejong_seoul_barplot_dodge_ggplotv2.png", bar_seoul_dodge, width = 6, height = 4.5, dpi = 72)

(bar_seoul_fill <- barplot.gg(vote.seoul.df, 
                              position = "fill", 
<<<<<<< HEAD
                              base_family = "HCR Dotum LVT", 
=======
                              base_family = "Malgun Gothic", 
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
                              ggtitle = "서울의 찬반", 
                              xlab = "계급", 
                              ylab = "집계(%)", 
                              fill.name = "찬반"))

<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

# ggsave("../pics/sejong_seoul_barplot_fill_ggplotv2.png", bar_seoul_fill, width = 6, height = 4.5, dpi = 72)
dump(c("barplot.gg","barplot.gg.stack", "barplot.gg.dodge", "barplot.gg.fill"), 
     file = "./barplot.gg.R")

mosaic plot

Data

vote.seoul.class.df <- as.data.frame(vote.seoul.class)
vote.seoul.class.sum <- tapply(vote.seoul.class.df$Freq, 
                               vote.seoul.class.df[, 2], 
                               sum)
vote.seoul.class.p.m <- prop.table(vote.seoul.class.sum)
vote.seoul.class.p <- prop.table(vote.seoul.class)
vote.seoul.class.p.2 <- prop.table(vote.seoul.class, 
                                   margin = 2)
vote.seoul.class.p.df <- as.data.frame(vote.seoul.class.p)
vote.seoul.class.p.df$width <- vote.seoul.class.p.m[match(vote.seoul.class.p.df$class, 
                                                          names(vote.seoul.class.p.m))]
vote.seoul.class.p.2.df <- as.data.frame(vote.seoul.class.p.2)
vote.seoul.class.p.df$height <- vote.seoul.class.p.2.df$Freq
vote.seoul.class.p.df$label.height <- unlist(tapply(vote.seoul.class.p.df$height,
                                                    vote.seoul.class.p.df$class, 
                                                    function(x) x / 2 + c(0, cumsum(head(x, -1)))))
vote.seoul.class.p.df$y.breaks <- unlist(tapply(vote.seoul.class.p.df$height,
                                                    vote.seoul.class.p.df$class, 
                                                    cumsum))
x.center <- vote.seoul.class.p.m / 2 + c(0, cumsum(head(vote.seoul.class.p.m, -1)))
# x.center <- (cumsum(vote.seoul.class.p.m) + 
#               c(0, head(cumsum(vote.seoul.class.p.m), -1)))/2
vote.seoul.class.p.df$center <- x.center[match(vote.seoul.class.p.df$class, 
                                               names(x.center))]

m1 <- ggplot(vote.seoul.class.p.df, 
             aes(x = center, 
                 y = height, 
                 width = width)) + 
  geom_bar(aes(fill = vote), 
           stat = "identity", 
           col = "white", 
           size = 2, 
           position = position_stack(reverse = TRUE))
m1

<<<<<<< HEAD

m2 <- m1 + 
  theme_bw(base_family = "HCR Dotum LVT")
m2

=======

m2 <- m1 + 
  theme_bw(base_family = "Malgun Gothic")
m2

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

m3 <- m2 + 
  geom_text(aes(x = center, 
                y = 1.05), 
            label = vote.seoul.class.p.df$class, 
<<<<<<< HEAD
            family = "HCR Dotum LVT")
m3

======= family = "Malgun Gothic") m3

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

m4 <- m3 + 
  geom_text(aes(x = center, 
                y = label.height), 
            label = format(vote.seoul.class.df$Freq, 
                           big.mark = ","), 
            position = position_identity())
m4

<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

x.breaks <- c(0, ifelse(cumsum(vote.seoul.class.p.m) < 0.1, 0.0, cumsum(vote.seoul.class.p.m)))
x.label <- format(x.breaks * 100, 
                  digits = 3, 
                  nsmall = 1)
y.breaks <- vote.seoul.class.p.df$y.breaks
delta <- (max(y.breaks) - min(y.breaks)) / 20
y.breaks.sort <- sort(y.breaks)
index <- which(diff(y.breaks.sort)  > delta)
y.breaks <- c(0, y.breaks.sort[c(index, length(y.breaks.sort))])
y.label <- format(y.breaks * 100,
                  digits = 2,
                  nsmall = 1)
m5 <- m4 + 
  scale_x_continuous(name = "계급(누적 %)", 
                     breaks = x.breaks, 
                     label = x.label) + 
  scale_y_continuous(name = "찬반(%)",
                     breaks = y.breaks,
                     label = y.label) + 
  scale_fill_manual(name = "찬반", 
                    values = rainbow(2)[2:1], 
                    labels = vote.seoul.df$vote, 
                    guide = guide_legend()) +
  ggtitle("서울의 찬반")
m5

<<<<<<< HEAD

=======

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

mosaic_gg

mosaic_gg <- function(tbl, base_family = "", 
                      ggtitle = "", 
                      xlab = "", 
                      ylab = "", 
                      fill.name = ""){
tbl.df <- as.data.frame(tbl)
N <- length(levels(tbl.df[, 1]))
tbl.sum <- tapply(tbl.df$Freq, tbl.df[, 2], sum)
tbl.p.m <- prop.table(tbl.sum)
tbl.p <- prop.table(tbl)
tbl.p.2 <- prop.table(tbl, margin = 2)
tbl.p.df <- as.data.frame(tbl.p)
tbl.p.df$width <- tbl.p.m[match(tbl.p.df[, 2], names(tbl.p.m))]
tbl.p.df$height <- as.data.frame(tbl.p.2)$Freq
tbl.p.df$label.height <- unlist(tapply(tbl.p.df$height, tbl.p.df[, 2], function(x) x / 2 + c(0, cumsum(head(x, -1)))))
tbl.p.df$y.breaks <- unlist(tapply(tbl.p.df$height, tbl.p.df[, 2], cumsum))
x.center <- tbl.p.m / 2 + c(0, cumsum(head(tbl.p.m, -1)))
# x.center <- (cumsum(tbl.p.m) + c(0, head(cumsum(tbl.p.m), -1)))/2
tbl.p.df$center <- x.center[match(tbl.p.df[, 2], names(x.center))]
m1 <- ggplot(tbl.p.df, aes(x = center, y = height, width = width)) + 
  geom_bar(aes(fill = tbl.df[1]), 
           stat = "identity", 
           col = "white", 
           size = 1, 
           position = position_stack(reverse = TRUE)) 
m2 <- m1 + 
  theme_bw(base_family = base_family)
m3 <- m2 + 
  geom_text(aes(x = center, y = 1.05), 
            label = tbl.p.df[, 2], 
            family = base_family)
m4 <- m3 + 
  geom_text(aes(x = center, y = label.height), 
            label = format(tbl.df$Freq, big.mark = ","), 
            position = position_identity())
x.breaks <- c(0, ifelse(cumsum(tbl.p.m) < 0.1, 0.0, cumsum(tbl.p.m)))
x.label <- format(x.breaks * 100, 
                  digits = 3, 
                  nsmall = 1)
y.breaks <- tbl.p.df$y.breaks
delta <- (max(y.breaks) - min(y.breaks)) / 20
y.breaks.sort <- sort(y.breaks)
diff(y.breaks.sort) < delta 
index <- which(diff(y.breaks.sort)  > delta)
y.breaks <- c(0, y.breaks.sort[c(index, length(y.breaks.sort))])
y.label <- format(y.breaks * 100,
                  digits = 2,
                  nsmall = 1)
m5 <- m4 + 
  scale_x_continuous(name = xlab, 
                     breaks = x.breaks, 
                     label = x.label) + 
  scale_y_continuous(name = ylab,
                     breaks = y.breaks,
                     label = y.label) + 
  scale_fill_manual(name = fill.name, 
                    values = rainbow(N)[N:1], 
                    labels = tbl.df[[1]], 
                    guide = guide_legend()) +
  ggtitle(ggtitle) +
  theme(plot.margin = unit(c(1, 2, 1, 1), "lines"))
return(m5)
}
dump(list = "mosaic_gg", file = "./mosaic_gg.R")

<<<<<<< HEAD

source("./mosaic_gg.R")
m5.seoul <- mosaic_gg(vote.seoul.class, 
                      base_family = "HCR Dotum LVT", 
=======
m5.seoul <- mosaic_gg(vote.seoul.class, 
                      base_family = "Malgun Gothic", 
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
                      ggtitle = "서울의 찬반", 
                      xlab = "계급", 
                      ylab = "찬반", 
                      fill.name = "찬반")
m5.seoul
<<<<<<< HEAD

ggsave("../pics/sejong_seoul_mosaic_ggplot.png", m5.seoul, width = 6, height = 4.5, dpi = 72)
=======

# ggsave("../pics/sejong_seoul_mosaic_ggplot.png", m5.seoul, width = 6, height = 4.5, dpi = 72)
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

지역별 관료들의 찬반(서울 제외)

kable(vote.region.bureaus[, -1], 
      caption = "지역별 관료들의 찬반(서울 제외)")

지역별 관료들의 찬반(서울 제외)

	경기	평안	황해	충청	강원	함길	경상	전라
찬성	29	6	17	35	5	3	55	42
반대	5	36	17	28	10	15	16	14

kable(format(prop.table(vote.region.bureaus[, -1], 
                        margin = 2)*100, 
             digits = 1, 
             nsmall = 1), 
      align = "r", 
      caption = "지역별 관료들의 찬반(서울 제외, %)")

지역별 관료들의 찬반(서울 제외, %)

	경기	평안	황해	충청	강원	함길	경상	전라
찬성	85.3	14.3	50.0	55.6	33.3	16.7	77.5	75.0
반대	14.7	85.7	50.0	44.4	66.7	83.3	22.5	25.0

막대그래프

vote.region.bureaus.df <- as.data.frame(vote.region.bureaus[, -1])
barplot.gg(vote.region.bureaus.df, 
           position = "stack", 
<<<<<<< HEAD
           base_family = "HCR Dotum LVT", 
=======
           base_family = "Malgun Gothic", 
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
           ggtitle = "지방 관료들의 찬반", 
           xlab = "지역", 
           ylab = "집계", 
           fill.name = "찬반")

<<<<<<< HEAD

ggsave("../pics/sejong_bureaus_barplot_stack_ggplot.png", width = 9, height = 4.5, dpi = 72)

=======

# ggsave("../pics/sejong_bureaus_barplot_stack_ggplot.png", width = 9, height = 4.5, dpi = 72)

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

mosaic plot

<<<<<<< HEAD

m5.bureaus <- mosaic_gg(vote.region.bureaus[, -1], 
                        base_family = "HCR Dotum LVT", 
=======
# source("./mosaic.gg.R")
m5.bureaus <- mosaic_gg(vote.region.bureaus[, -1], 
                        base_family = "Malgun Gothic", 
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
                        ggtitle = "지방관료의 찬반", 
                        xlab = "지역별 관료들의 비중(누적 %)", 
                        ylab = "찬반(%)", 
                        fill.name = "찬반")
m5.bureaus
<<<<<<< HEAD

ggsave("../pics/sejong_bureaus_mosaic_ggplot.png", m5.bureaus, width = 9, height = 4.5, dpi = 72)
=======

# ggsave("../pics/sejong_bureaus_mosaic_ggplot.png", m5.bureaus, width = 9, height = 4.5, dpi = 72)
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

품관촌민들의 지역별 찬반

막대그래프

vote.region.commons.df <- as.data.frame(vote.region.commons)
barplot.gg(vote.region.commons.df, 
           position = "stack", 
<<<<<<< HEAD
           base_family = "HCR Dotum LVT", 
=======
           base_family = "Malgun Gothic", 
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
           ggtitle = "품관촌민의 찬반", 
           xlab = "지역", 
           ylab = "집계", 
           fill.name = "찬반")

<<<<<<< HEAD

ggsave("../pics/sejong_commons_barplot_stack_ggplot.png", width = 9, height = 4.5, dpi = 72)

=======

# ggsave("../pics/sejong_commons_barplot_stack_ggplot.png", width = 9, height = 4.5, dpi = 72)

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

mosaic plot

<<<<<<< HEAD

m5.commons <- mosaic_gg(vote.region.commons, 
                        base_family = "HCR Dotum LVT", 
=======
# source("./mosaic.gg.R")
m5.commons <- mosaic_gg(vote.region.commons, 
                        base_family = "Malgun Gothic", 
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
                        ggtitle = "품관촌민의 찬반", 
                        xlab = "지역별 품관촌민의 비중(누적%)", 
                        ylab = "찬반(%)", 
                        fill.name = "찬반")
m5.commons
<<<<<<< HEAD

ggsave("../pics/sejong_commons_mosaic_ggplot.png", m5.commons, width = 12, height = 4.5, dpi = 72)
=======

# ggsave("../pics/sejong_commons_mosaic_ggplot.png", m5.commons, width = 12, height = 4.5, dpi = 72)
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

충청도

관료들의 찬반과 품관촌민의 찬반이 다른 곳.

kable(format(vote.chung.class, 
             big.mark = ","), 
      caption = "충청도 사례", 
      align = "r")

충청도 사례

	대신 등	수령	품관촌민
찬성	0	35	6,982
반대	2	26	14,013

kable(format(prop.table(vote.chung.class, 
                        margin = 2)*100, 
             digits = 3, nsmall = 1), 
      caption = "충청도 사례", 
      align = "r")

충청도 사례

	대신 등	수령	품관촌민
찬성	0.0	57.4	33.3
반대	100.0	42.6	66.7

막대그래프

vote.chung.class.df <- as.data.frame(vote.chung.class)
barplot.gg(vote.chung.class.df, 
           position = "fill", 
<<<<<<< HEAD
           base_family = "HCR Dotum LVT", 
=======
           base_family = "Malgun Gothic", 
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
           ggtitle = "충청도의 찬반", 
           xlab = "계급", 
           ylab = "집계(%)", 
           fill.name = "찬반")

<<<<<<< HEAD

ggsave("../pics/sejong_chungcheong_barplot_fill_ggplot.png", width = 6, height = 4.5, dpi = 72)

=======

# ggsave("../pics/sejong_chungcheong_barplot_fill_ggplot.png", width = 6, height = 4.5, dpi = 72)

>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

mosaic plot

<<<<<<< HEAD

m5.chungcheong <- mosaic_gg(vote.chung.class[, -1], 
                            base_family = "HCR Dotum LVT", 
=======
# source("./mosaic.gg.R")
m5.chungcheong <- mosaic_gg(vote.chung.class[, -1], 
                            base_family = "Malgun Gothic", 
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86
                            ggtitle = "충청도의 찬반", 
                            xlab = "계급(누적%)", 
                            ylab = "찬반(%)", 
                            fill.name = "찬반")
m5.chungcheong
<<<<<<< HEAD

ggsave("../pics/sejong_chungcheong_mosaic_ggplot.png", m5.chungcheong, width = 12, height = 4.5, dpi = 72)
=======

# ggsave("../pics/sejong_chungcheong_mosaic_ggplot.png", m5.chungcheong, width = 12, height = 4.5, dpi = 72)
>>>>>>> f5332224c7063b88292b34e1d129b95716d0fa86

자료 정리

save.image(file = "./sejong_ggplot.RData")